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Abstract 

We evaluate numerically several superconducting correlation functions in a 
generalized t — J model derived for hole-doped Cu02 planes. The model 
includes a three-site term t" similar to that obtained in the large U limit 
of the Hubbard model but of opposite sign for realistic 0-0 hopping. For 
realistic parameters we obtain strong evidence of superconductivity of pre- 
dominantly d x 2_ y 2 character. The ground state has a large overlap with a 
very simple resonating-valence-bond wave function with off-diagonal long- 
range order. This function reproduces the main features of the magnetic 
and superconducting correlation functions. 
PACS numbers: 74.20.Mn, 57.10.Jm, 75.40.Mg 
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In spite of a considerable effort, no indisputable evidence of superconductivity has been 
found in strongly correlated models for the cuprates for realistic parameters. The numerical 
methods are so far the most reliable to treat strong correlations, and clear indications of 
d x 2_ y 2 superconductivity have been found in the t — J model for a 4 x 4 cluster, but only for 
values of J which are an order of magnitude larger than the realistic values and for too large 
doping x = 0.5 jy. Variational Monte Carlo results in the large U limit of the Hubbard 
model suggest that it might be necessary to solve larger clusters to obtain superconductivity 
for realistic parameters, and favor mixed s — d superconductivity ||. On the other hand, 
Monte Carlo results in the three-band Hubbard model |||4[] with 0-0 hopping t pp = 
favor s-wave superconductivity, although the dependence of the superconducting correlation 
functions with distance has not been investigated 0. In this Letter we present numerical 
evidence of superconductivity and at the same time a resonating-valence-bond (RVB) state 
in a realistic model. 

After the original derivation of the t — J model [0], it became clear that this model 
should be supplemented by other terms to accurately represent the physics of the three- 
band Hubbard model [PHT2|, and experimental data |T2|-|T4|]. Including the most important 
terms, the Hamiltonian can be written in the form fll| |: 
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H acts on a square lattice on which double occupancy is forbidden (c\ a rii = 0). The vectors 
5(7) connect a site with each of its four nearest neighbors (next-nearest neighbors). While 
the sign of t can be absorbed in a change of phase of half the creation operators, the signs 
of t' and t" are relevant and depend on the form chosen for the restricted Hilbert space and 
the precise meaning of cJ CT . Here this operator creates a particle at the vacant site i [ITS . 



A recent systematic analytical study of the mapping of the three-band Hubbard model 



to H shows that for 0-0 hopping t pp > 0.3eV, the mapping using localized non-orthogonal 
Zhang-Rice singlets is more accurate than the one which uses orthogonal O Wannier func- 
tions ||12|| . Using the former mapping, t pp gives rise to a positive contribution to t" due to the 
non-orthogonality of the states. As shown in Fig. 1, numerical fitting of the parameters of 
H agrees very well with the analytical calculations and both show that t' and t" change sign 
for increasing t pp . Constrained density-functional calculations predict t pp =0.65eV || and 
this results in positive values of t' and t". A positive t' [15|] agrees with calculations of other 



authors [§,0], with the form of the Fermi surface JT4J, and with calculations of different 



normal state properties |[L3| . Also, the fact that i! > for hole-doped systems while t' < 



in electron doped systems |9|,[r|,|14]] allows to understand the different relative stability of 



the Neel order upon doping in both types of systems [pL2|,p~4|| - Thus, the parameters taken 
here: J = 0.4, t' ~ 0.2, and t" ~ 0.1 in units of t = 1, are close to those which best describe 
the hole-doped cuprates. 

The effect expected from the correlated hopping t" becomes clear writing the correspond- 
ing term in the form of a hopping of nearest-neighbor singlets: 

H t „ = 2t" £ b\ 6 ,b id ; b\ 5 = -U4 +(5 T c i| - cU4 T ) • (2) 

Thus, H t n is expected to favor a Bose condensate of singlets with d x 2_ y 2 (s) symmetry for 
t" > 0(t" < 0), similar to the RVB state constructed by Anderson 0. This agrees with 
mean-field calculations 0], and the above mentioned results for the Hubbard model (for 



which t" = — J I '4) and the three-band Hubbard model with t pp — || (implying also t" < 
in the effective model H). Instead J at most tends to stabilize static singlets, while t and 
part of the terms proportional to t', break the singlets. Since c?-wave pairing is expected 
to arise from antiferromagnetic spin fluctuations at small doping flI||PT||, one expects that 
a positive t" would interfere constructively with any superconducting mechanism already 
present in H — H t " . 

We have used the Lanczos method to obtain the ground state of H for a 4 x 4 cluster with 
periodic boundary conditions in each subspace defined by the z projection of the total spin 



and the irreducible representation of the space group. In Fig. 2 we show the susceptibility 
xtup = (T,i<j68> /(<*)/(<*') < 9 I b\s b js> I 9 >)/L where | g > is the ground state, L the 
size of the system, f(5) — 1 if 5 || x and /(<5) = — 1 if 5 || y as a function of 
The abrupt increase in xtup ls apparent. After this increase, x d sup reaches values similar 
or larger than the maximum one obtained previously in the t — J model (xtup ~ 2.5 for 
J = 3,x = 0.5, t' = t" = 0) Q. Also, the qualitative behavior of the superconducting 
correlation functions (shown below) changes and they saturate rapidly with distance for 
x = 0.25 but not for x > 0.5. As will be described in detail elsewhere, except for very 
particular values of the parameters, xtup has its maximum value at a doping level x ~ 0.25 
in agreement with experiment. For x = 0.75, x d SU p ~ 10~ 2 f° r ^ = 0-2- F° r all parameters 
considered here the system is far from the region of phase separation 0, which is reduced 
by t". 

Our results suggest that as t" is increased, there is a continuous, although abrupt tran- 
sition from a disordered spin liquid to a RVB state with off-diagonal long range order 
(ODLRO). To support this conclusion we consider the following simple generalization to 
x ^ and ODLRO of Sutherland's RVB wave function @: 

\RVB0>=P N Y[(l + Y / g(5)b\ s )\0>, (3) 

jeA S 

where Pjsr is the projector over a definite even number of particles N, jeA indicates that the 
product over sites is restricted to one sublattice g(5) = 1 for each "horizontal" singlet (<5 || x) 
and g(5) = i for each "vertical " singlet (8 || y). Eq. (3) is clearly not an optimum wave 
function, since unlike Anderson's one it contains no wave-vector dependence. In spite of 
this we find that from its real part and the first few powers of H applied to it, an excellent 
variational function for t" > 0.12 can be constructed. The real (imaginary) part of | RVB0 > 
has an even (odd) number of vertical singlets. For N = L, Im | RVB0 >=0. It is easy to 
check that if the number of singlets N/2 is a multiple of four, Re \ RVB0 > transforms like 
the representation r x (invariant) of the space group, while Im | RVB0 > transforms like 
T3 (x 2 — y 2 ). Instead for N/2 even but not multiple of four, the real (imaginary) part of 
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I RVBO > transforms like ^ (Ti) under the symmetry operations of the space group. In 
the 4 x 4 cluster we find that for N = 12 (x = 0.25) and N = 8 (x = 0.5) for t" > 0.1, the 
ground state | g > belongs to the same irreducible representation as Re \ RVBO >. The 
simplest function which contains configurations with even and odd number of vertical bonds 
with the correct symmetry, and takes advantage of H t /> has the form: 

| RVB1 >= F(l - jH t „)Re \ RVBO > , (4) 

where F is a normalization factor. We have determined a maximizing the square of the 
overlap S 2 =\< g \ RVB1 >\ 2 . This quantity as a function of t" is shown in Fig. 3. Taking 
into account the simplicity of Eq. (4), that there is only one free parameter, and that the 
size of the Hilbert space is huge, the magnitude of S 2 after the transition is noticeable. 

| RVB1 > can be improved considerably if a variational function is constructed with 
H n | RVB1 > {n = to 3). The generalized RVB state thus obtained correctly takes into 
account the effects of t and t' and practically coincides with | g >. However, we will show 
that | RVB1 > reproduces already, at least qualitatively, the main numerical results. xtu P 
obtained from | RVB1 > is compared in Fig. 2. Note that the expectation value of b\ s biS' 
in | RVB1 > for 5 _L 5' vanishes for a = 0, implying equal amounts of s— and <i-wave 
superconductivity. For t" > 0.12 the <i-wave component of | RVB1 > dominates but there 
is always some s component. 

In Fig. 4 we show how after the transition with increasing t", a peak in (7r, 7r) develops in 
the magnetic structure factor, which is in excellent agreement with the result of | RVB1 >. 
This peak is probably related with the strong peak in the staggered spin susceptibility x( n , n ) 
assumed in the phenomenological spin-fluctuation theories |i~7| . For t" = 0, S(tt, it) slightly 
decreases with increasing t'. The absence of a structure in S(ir/2,7r/2) for | RVB1 >, is 
due to the lack of magnetic correlations beyond nearest neighbors in this function, and is 
restored if | RVB1 > is improved including the effect of t. 

In Fig. 5 we show the distance dependence of the d-w&ve superconducting correlation 
functions c(m) = J2i8&' f(d)f(8') < 9 I bj s b i+m s> \ g > and the corresponding expectation 
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value in | RVB1 >. One can see that c(m) saturates at large distances, suggesting ODLRO 
for the infinite system. The significant overlap of | g > with a function which has ODLRO 
in the thermodynamic limit is another important indication. Instead, the corresponding 
s-wave correlation functions decay steadily with distance. 

In summary, we have found strong evidence of <i-wave superconductivity in the model 
Eq. (1) for reasonable t" > 0.12 and experimentally relevant doping x = 0.25, and parameters 
J = 0.4, t' ~ 0.2 in units of t. The possibility of some amount of s-wave superconductiv- 
ity remains open. The model and parameters are close to the most realistic ones for the 
cuprates. The superconductivity is rather insensitive to if and preliminary results suggest 
that it persists if a repulsion between nearest-neighbor Zhang-Rice singlets smaller than J 
is included in the model. The ground state of the system for the above mentioned parame- 
ters is a generalized superconducting RVB state || . A similar conclusion has been recently 
obtained for single-rung t — J ladders |19j . As noted by Anderson, the RVB state is favored 
by electron-phonon interaction, particularly if the system is near an instability against a Cu 
dimerization mode ||. Theoretical calculations have suggested that this might well be the 
,PTJ. Also a RVB state allows for a qualitative explanation of the phase diagram 



case 



of the high-Tc systems at mean-field level pifl , and is consistent with a strongly enhanced 
x(7r, 7r) [17| and with the spin gap observed in YBaCuO [23|. To our knowledge this is the 
first time that large superconducting correlations at distances of a few lattice parameters 
and the physics of the RVB is unambiguously obtained numerically in the ground state of a 
model for hole-doped Cu02 planes for parameters close to the optimum ones. 

We would like to thank M. Balina for providing us with useful numerical subroutines 
he developed, and to M.E. Simon, M. Balina, E. Gagliano, B. Alascio and E. Dagotto for 
useful discussions. One of us (CDB) is supported by the Consejo Nacional de Investigaciones 
Cientfficas y Tecnicas (CONICET), Argentina. (AAA) is partially supported by CONICET. 



6 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 
[7 



[9 

[10 

[11 
[12 

[13 



[14 



REFERENCES 

E. Dagotto and J. Riera, Phys. Rev. Lett. 70, 682 (1993). 

G.J. Chen, R. Joynt, F.C. Zhang, and C. Gros, Phys. Rev. B 42, 2662 (1990); Q.P Li, 
B.E.C. Koltenbah, and R. Joynt, Phys. Rev. B 48, 437 (1993). 

V.J. Emery, Phys. Rev. Lett. 58, 2794 (1987). 

P.B. Littlewood, CM. Varma, and E. Abrahams, Phys. Rev. Lett. 63, 2602 (1989) and 
references therein. 

G. Dopf, A. Muramatsu, and W. Hanke, Phys. Rev. Lett. 68, 353 (1992). 
PW. Anderson, Science 235, 1196 (1987). 

F. C. Zhang and T.M. Rice, Phys. Rev. B 37, 3759 (1988). 

R.J. Gooding and V. Elser, Phys. Rev. B 41, 2557 (1990); F.C. Zhang and T.M. Rice, 
ibid 2560. 

M.S. Hybertsen, E.B. Stechel, M. Schliiter, and D.R. Jennison, Phys. Rev. B 41, 11068 
(1990). 

H. Q. Ding, G.H. Lang, and W.A. Goodard III, Phys. Rev. B 46, 14317 (1992). 

CD. Batista and A. A. Aligia, Phys. Rev. B 48, 4212 (1993); 49, 6436 (E) (1994). 

A.A. Aligia, M.E. Simon, and CD. Batista, Phys. Rev. B. 49, 13061 (1994) 

H. Chi and A.D.S. Nagi, Phys. Rev. B 46, 421 (1992); M.J. Lercher and J.M. Wheatly, 
Physica C 215, 145 (1993); R.J. Gooding, K.J.E. Vos, and PW. Leung, Phys. Rev. B 
49, 4119 (1994). 

T.Tohyama and S. Maekawa, Phys. Rev. B 49, 3596 (1994). 



[15] If the Zhang- Rice singlets were represented by doubly occupied sites as in Ref. 0, or 
if the c\ a would destroy particles creating vacant sites, all tfs would change sign. In 



most studies of t — J- like models the particles are called "electrons" and the vacant 
sites "holes", although for hole-doped systems they represent respectively holes and 
Zhang-Rice singlets of the original multiband Hubbard model. 

[16] G. Kotliar, Phys. Rev. B 37, 3664 (1988). 

[17] N.E. Bickers D.J. Scalapino, and S.R. White, Phys. Rev. Lett. 62, 961 (1989); P. Mon- 
thoux, A.V. Balatsky, and D. Pines, Phys. Rev. Lett. 67, 3448 (1991); P. Monthoux 
and D. Pines, Phys. Rev. Lett. 69, 961 (1992). 

[18] N. Read and B. Chakraborty, Phys. Rev. B. 40, 7133 (1989); and references therein. 

[19] S. Gopalan, T. M. Rice and M. Sigrist, Phys. Rev. B 49, 8901 (1988); M. Sigrist, T. M. 
Rice and F. C. Zhang, Phys. Rev. B (to be published). 

[20] S. Barisic, I. Batistic, and J. Friedel, Europhys. Lett. 3, 1231 (1987). 

[21] A.A. Aligia, Phys. Rev. B 39, 6700 (1989). 

[22] PA. Lee and N. Nagaosa, Phys. Rev. B 46, 5621 (1992). 

[23] J. Rossat-Mignod et al, Physica C 185-189, 86 (1991); Physica B 186-188, 1 (1993). 



8 



FIGURE CAPTIONS 



Fig.l: Parameters of H determined from the three-band Hubbard model as a function of 
the 0-0 hopping t pp . The remaining parameters were taken from Ref. |J and the Cu-0 
hopping t p d was taken as the unit of energy. Full line: numerical fit of the energy levels 
of a CU4O8 cluster. Dashed line: analytical results from a mapping using localized non- 
orthogonal Zhang-Rice singlets. 

Fig. 2: (i-wave susceptibility as a function of t" for t — 1, J — 0.4 and doping x = 0.25. Full 
line: t' = 0.2 . Dashed line: t' = 0. Dotted line: result using | RVB1 > (Eq. (4)) with 
t' = 0.2. 

Fig. 3: Square of the overlap between the ground state and | RVB1 > as a function of t" 
for t' = 0.2. Other parameters as in Fig. 2. 

Fig. 4: Magnetic structure factor < 9 I Sf 5Je ?q ^ Ri_R; '^ | g > /L as a function of 
wave vector. Dashed line: t' = t" = 0. Full line: t' = t" = 0.2. Dotted line: result using 
I RVB1 > with t' = t" = 0.2. Other parameters as in Fig. 2. 

Fig. 5: rf-wave superconducting correlation functions as functions of the distance. Param- 
eters are the same as in Fig. 4. 
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